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Hysteresis curves for finite arrays of JV x JV ferromagnetic nano dots subject to the dipole-dipole 
interaction are investigated for JV = 2 ... 13. Spin arrangements up to JV = 6 are presented, which 
indicate the onset of bulk-like behavior associated with odd (TV = 5) and even (JV = 6) systems. 
The effect of field misalignment on the hysteresis loops is also studied for JV = 3 • ■ • 6. The area An 
of the hysteresis loop is studied as a function of JV. We find that An — Aoo approximately scales as 
JV~2 for JV odd and as JV~ 2 for JV even. 
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INTRODUCTION 

A ferromagnetic particle goes into a monodomain state 
if its size D is below a critical value D c = 10 ~ 100 nm. 
This is due to the competition between the exchange and 
dipolar energies. Therefore, a nanoparticle in a mon- 
odomain state may be viewed as a giant magnetic dipole 
with magnetic moment of thousands of Bohr magnetons. 
For an JV x JV array of well-separated nanoparticles the 
exchange energy is usually negligible in comparison with 
the dipolar and anisotropy and Zeeman energies. The 
study of such systems is of increasing importance be- 
cause of their technological applications in data storage 
devices and magnetic field sensors. As the technology of 
these devices moves towards higher densities of stored in- 
formation, it requires smaller particles of magnetic media 
0, 0, for which finite size effects become relevant. In 
finite arrays of such large dipole moment particles, the 
dipolar field of the array becomes comparable with the 
bulk anisotropy field. Dipolar effects in such systems af- 
fect the static and dynamics properties of the array; and 
thus must be taken into account. 

Recently, Camley and Stamps in Q, Q, || investi- 
gated the dynamics and magnetization processes of a fi- 
nite planar array of JV x JV ferromagnetic nano dots, for 
JV = 3, 4, 5, 6. The nano-dots were taken to interact only 
via the dipole-dipole interaction, and they were subject 
to an external field applied either along one side of the 
array or along its diagonal. They found rather compli- 
cated hysteresis loops with the magnetization reversal 
controlled by the shape anisotropy induced by the array 
itself. We have considered the same model, and have ex- 
tended their results, for JV — 2 • • • 13. Our results for 
JV = 3 qualitatively agree with those of [Si]. We find 
that the behavior of these systems is surprisingly com- 
plex, both for small and for larger values of JV, and we 
have studied the approach of these systems to JV — ► oo 
behavior. 

In the present work each dot is taken to have a ra- 
dius Rd, thickness d and a single degree of freedom cor- 
responding to the orientation of a magnet of saturation 
magnetization Mq. We consider only the case of zero 



temperature. The dots are arranged on a square lattice 
with lattice spacing a > 2Rd, and the dots interact only 
via the dipole-dipole interaction. The equation of motion 
for the magnetic moment of each dot is governed by the 
Landau-Lifshitz-Gilbert equation (LLG), 6] which reads 
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where 7 is the gyromagnetic ratio, a is the damping 
coefficient, M is the magnetic moment of the dot and 
M s = |M| is the saturation magnetization, H e /r is the 
average effective magnetic field acting on the dot. The 
average effective magnetic field acting on the i-th dot is 
due to the applied external field, the dipolar fields, and 
the anisotropy field 

H* e// = H cos 9x + H sin 9y - W dip + ^K^z. (2) 

Here the dipole field acting on the i-th dot due to all 
other dots in the array is given by 
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The choice of anisotropy field is determined by the 
shape of the dot, which in our problem is directed along 
the symmetry axis of the (cylindrical) dots. We divide 
both sides of Eq.JTJ by (7M5) and define a dimensionless 
time variable r = jM s t. The LLG in these reduced units 
becomes 



— = m x h e// - -m x (m x h e// ) (4) 
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We employ a system of 



where m - M<> _ ^ eJJ - M> 
units where magnetic fields are given in units of M s and 
distances are given in units of the array's lattice spacing 
a. The strength of the dipole field is characterized by 
hdip = ^7p~ 1 which is the ratio of the dot volume to the 
volume of a cubic cell with side a. For all arrays studied 
in this work, we take hdip = 0.5. 
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This article is organized as follows. Section II presents 
the numerical techniques and method of time integration. 
Section III presents an extensive discussion of magneti- 
zation processes and hysteresis for arrays of A x TV nano 
dots (A — 2, 3, • • • , 13) when the external field is applied 
along one side of the array. Section IV considers the 
effects on the hysteresis loop of misalignment of the ex- 
ternal field. Section V considers the relationship between 
the area of the hysteresis loop and N. A brief summary 
is given in section VI. 

NUMERICAL TECHNIQUES 

We have employed two different approaches to study 
the magnetization processes of our A x A arrays of nano- 
dots. The first approach uses the second rank Runge- 
Kutta (RK) algorithm with fixed time step to integrate 
the LLG equation in reduced units. The second approach 
uses the "greedy algorithm" to find the stable final state 
of the array. In both approaches the calculations were 
done using FORTRAN and Mathematica languages 0. 
The two approaches yielded similar results. 

In the RK approach, the integration is done using a 
fixed time step At = 5 x I0~ 3 and a damping coefficient 
— = 0.6, with an initial state in which the magnetic 
moment of the dot is randomly generated. The time 
integration proceeds until a stable final state is reached. 
The stability of the final state is checked by changes in the 
total energy of the system. Iterations are stopped when 
the difference between the total energy of the system from 
the (n+ l)-th iteration and that of the n-th iteration is of 
the order of AE n = 10~ 5 in units of M^a 3 . Our solutions 
converged after almost 10 3 iterations. 

The greedy algorithm approach assumes the dots to 
be aligned along the direction of their total local field. 
In the initial state each dot is chosen to point randomly. 
Next, all components of the total local field are calcu- 
lated for each dot. Finally each component of the field 
is normalized by the magnitude of the total local field to 
yield the new dot magnetization orientation. The final 
state of the n-th iteration is used as an initial one for the 
next iteration until the solution is converged. The con- 
vergence of the final state is checked in manner similar 
to that used in the first approach. Note that the greedy 
algorithm can be obtained from the LLG equation in the 
limit where the damping coefficient is very large. We find 
that the LLG approach converges faster than the greedy 
algorithm. This is probably related to a phenomenon 
known in RLC circuits, where a critically-damped cir- 
cuit approaches equilibrium faster than a circuit with a 
larger amount of damping. 

Calculation of the dipolar field at the z-th dot is the 
most computationally time-consuming aspect of both ap- 
proaches since it requires summation of the dipole fields 
from all other dots in the array. However, due to the 



relatively small sizes of our dot arrays this calculation is 
performed rather quickly. 

HYSTERESIS LOOP AND EXTERNAL FIELD 
ORIENTATION EFFECT 

Hysteresis loops M(H) for A x TV arrays of nano dots 
subject to an external magnetic field applied along one 
side have been calculated. Initially, a strong external field 
Hq is applied to the array until saturation. Then the field 
is decreased to —Hq, and then is increased again up to 
Hq. We have taken Hq = 2M S , and a fixed field-step 
of AH = 2 x 10 _3 M S is used to simulate the sweeping 
process. For each value of the external field the system 
was iterated until a stable final state is reached. As seen 
in Fig. ^ where M(H) is plotted (both in units of M s ), 
the odd and even TV arrays have somewhat different be- 
haviors, especially for small A. One aspect of this is that 
the odd-A systems display magnetization jumps as the 
field changes. The odd-TV and even-TV behavior becomes 
similar for larger values of A, something we study in a 
later section. 

In Fig.QJ the angle of the field to one of the sides (the 
y-axis) is taken to be 8 = 0. Experimentally, however, 
field misalignment is almost inevitable, so that we have 
also studied this phenomenon. 

Fig. shows results for N = 3 ... 6 and angles 8 = 5°, 
30°, and 45°. (We present only some of the more repre- 
sentative results; angles between and 45° were studied 
in 5° increments.) 

Comparison of Fig. with Fig.|3for 8 = 5° shows that 
a small misalignment of the applied field can change the 
hysteresis loop drastically. 

For N = 3, Fig|21 shows that the central part of the 
hysteresis loop shrinks as 8 increases. For 8 = 45°, the 
central part almost disappears completely, and new small 
loops start to develop away from the center of the hys- 
teresis loop. Our results for A = 3 agree with those given 
m 

3], 

which studied the cases 8 = 0° and 45°. 

For TV = 4, at 8 = 0° there is no central loop, but there 
is a prominent loop at finite field. On the other hand, 
at 8 = 5°, there is a central loop, and the finite field 
structure becomes rather complex. Further increase of 8 
leads to a filling out and connecting of various subloops. 
Also, note the appearance of jumps for non-zero 8. 

For A = 5, a small field misalignment has an enormous 
effect, at 8 = 5° shrinking the loop to a relatively small 
central region. As 8 increases, the central loop grows, but 
the loop for 8 = 45° pinches off to yield three subloops, 
as for A = 3. 

For TV = 6, again a small field misalignment has an 
enormous effect, at 8 = 5° shrinking the loop to a rel- 
atively small central region. As 8 increases, the central 
loop grows, but in contrast to A = 5, the loop for 8 = 45° 
does not pinch off, and closely resembles the loop for 
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FIG. 1: Hysteresis loops M(H), both in units of M a , for 
weakly coupled arrays of N x N ferromagnetic nano dots. 
The external field is applied along the y-axis, which coincides 
with one side of the array. The top two rows are for N even 
and the lower two are for N odd. 

N = 4. 

These different types of behavior indicates that these 
are complex systems, for which it is difficult to make 
generalizations. 

HYSTERESIS AND EVEN-ODD SIGNATURE IN 
FINITE ARRAY OF NANO-DOTS 

In the absence of an external magnetic field the array 
of N x N nano-dots favors antiferromagnetic ordering, 
thus minimizing its magnetostatic energy. A large exter- 
nal magnetic field applied to the array tends to orient the 
magnetic moments along the field, thus minimizing the 
Zeeman field energy. However, the spins at the array cor- 
ners tip by a small angle, as shown in Fig. El-Fig. El form- 
ing a two-dimensional "flower" state [8(, |9|- The flower 
state persists until the applied field falls to Hq = M s . For 
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FIG. 2: Hysteresis loops M(H), both in units of M 3 , for 
N = 3 • • • 6 arrays for the field at 8 = 5°, 30° and 45° going 
from left to right. 

lower values of the applied field, the competition between 
the dipole-dipole interaction and the Zeeman energy be- 
comes significant, and changes the array ordering. 

For N = 2, Fig. shows that below Hq = M s and 
for Hq ^ the spins form a snake-like domain structure 
winding clockwise or counterclockwise. At Hq = 0, the 
array has zero net per-dot magnetization, due to a vortex- 
like structure that persists for — 3M S < Hq < 0.3M S . 

The N = 3 array was analysed by Camley and Stamps 
in Q . This array also shows snake- like arrangements be- 
low Hq = M s . At zero field the final state of the array 
shows what we call a " barrel" state in which the spins at 
the left and right columns are oriented opposite to the 
central column with a slight tipping of the corner spins, 
as shown in Fig. 0] This agrees with Ref. except for the 
tipping of the corner spins. We have made a small angle 
expansion of the energy for the spins being nearly aligned, 
and indeed we find that the tipped corner spin state is 
more energy favorable. The tipping angle was determined 



H=2.0 



H=0.8 



H=0.5 




H=0.3 



H=0.0 



H=-0.3 



H=-0.5 



H=-0.8 



H=-2.0 



FIG. 3: Spin arrangements for an array of 2 x 2 ferromagnetic 
nano dots in external magnetic field. 
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FIG. 5: Spin arrangements for an array of 4 x 4 ferromagnetic 
nano dots in external magnetic field. 



H=2.0 



t ! 
; t 



i ! 

i ! 

1 1 

1 t 

\ \ 

H=0.3 



\ 


^ X 


/ \ / 


/ / 


/ \ / 


/ / 


\ /\ 
/^-- 


v \ 
W 


H= 


=-0.5 


/ 


y v 


\ \ / 


\ / 


1 I I 

I I \ 


I 1 

/ \ 


\ x^ 





H=0.8 


; 


\\\ 1 
/ / / / 


! 

\ 

/ 


tilt 
\ \ \ \ 
// / | 


H=0.0 


\ 




/ 


\ / \ f 


1 
\ 


I 1 I I 

/ \ / \ 


/ 




H=-0.8 




/ / / / 
WW 


I 


I 1 I J 


1 


/ / / I 


\ 


WW 



H=0.5 



t / ^ / / 

H=-0.3 



1 I I 


'H 











H=-2.0 



\ 1 


i 1 


I I 
1 I 


I I 
I 1 


I I 


I 1 


I I 


i \ 



FIG. 4: Spin arrangements for an array of 3 x 3 ferromagnetic FIG. 6: Spin arrangements for an array of 5 x 5 ferromagnetic 
nano dots in an external magnetic field. nano dots in external magnetic field. 



by both iteration and analytical calculations, and agree- 
ment is found between the two values. The numerical 
value of the corner spin tipping angle, found both by it- 
eration and small tipping angle analysis, is |cv| = 9.1115°, 
with numerical error of order 10~ 5 . The spin snapshots 
in Fig. 21 and hysteresis loop analysis show that the bar- 
rel state switches to an inverted barrel state when the 



applied field changes sign. 

Fig[S] shows that the N = 4 array also features snake- 
like arrangements of the spins, for intermediate values of 
the applied field. However, in zero field the total per- 
dot magnetization is zero, which can be attributed to 
the formation of a vortex in the array's central 2x2 
block. The magnetic moments of the rest of the dots 
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FIG. 7: Spin arrangements for an array of 6 x 6 ferromagnetic 
nano dots in external magnetic field. 



in the array form a ring that surrounds the vortex with 
opposite circulation. This state is stable for applied fields 
H satisfying -0.2M S < H < 0.2M S . 

For all TV, the flower state appear at high fields (here, 
| lio | = 2.0M S ). The hysteresis loops shown in Fig. 
show a subtle difference in shape between arrays with 
odd TV and arrays with even TV. For odd TV the loops 
show well-defined jumps whereas for even TV this behav- 
ior is absent. This behavior is due to unpaired spins 
with uncompensated dipole fields. The jumps become 
less apparent as TV grows, and will eventually disappear 
for large TV, where the distinction between even and odd 
TV becomes unimportant. 

Note that, when an array was placed in zero exter- 
nal field and given random initial conditions, the solu- 
tions converged to the same states as obtained in the 
hysteresis-cycle calculations, up to the degeneracy of the 
system. Thus, for TV — 2,4 there are two degenerate 
metastable states of opposite chirality (winding) with 
zero net magnetization in zero field, each of which has 
a four-fold rotational symmetry. For TV — 6 there are 
two degenerate states of opposite chirality, with non- 
vanishing net magnetization, each of these state has no 
apparent rotational symmetry. For TV = 3 there are two 
degenerate barrel states, with no rotational symmetry, 
and TV = 5 is similar to TV = 6. The hysteresis loop area 
An will produce further evidence that large system be- 
havior commences with TV = 5 and TV = 6. Metastable 
states with vanishing net magnetization may appear for 
arrays with even TV. However, our simulations showed 
that these states could appear only for TV — 2, 4. For 



arrays with odd TV the unpaired dipoles prevent the for- 
mation of such states. 



HYSTERESIS LOOP AREA A N VS. PARTICLE 
NUMBER TV 



Although the area of the hysteresis loop An tends to 
zero for the TV = 2 array, it clearly is nonzero for all other 
arrays. We present the hysteresis loop areas in Fig. [SI as 
circles and squares. From Fig. |H1 we see that the area 
of the hysteresis loop decreases with increasing TV except 
for TV = 3 for TV odd and TV = 2,4 for TV even. The 
TV — 5 and TV = 6 arrays, the first to show something 
like bulk behavior, have maximum An for odd and even 
respectively; their spin arrangements are given in Fig. 
and Fig. [7| We have fitted our data to the asymptotic 
form 



.4 



N 



C_ 

NP 



(5) 



where A^, C and p are constants to be determined. If 
larger values of TV had been computationally feasible, we 
would have considered only large values of TV for the fit. 
The fit is not as good when TV = 5 is included, so we do 
not show this case. In practice, for TV odd the data are fit 
starting from TV = 7 and for TV even from TV = 6. Both 
fits are shown as solid lines in Fig. El where the values of 
Aoo, C, and p are given in the table. 




FIG. 8: The area of the hysteresis loop as a function of the 
number of particles TV. 

The associated values of x 2 are both less than 10~ 5 . 
For odd TV, An varies approximately as An ~ TV" 2 
whereas for even TV it varies approximately as TV" 2 . We 
attribute no fundamental significance to these values. 
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TABLE I: Fitting parameters for data given in Fig. |H| using 
Eq.EI 



N 


^4-00 


C 


P 


even 


0.278 


6.31 ±0.42 


1.95 ±0.03 


odd 


0.278 


6.92 ±0.86 


1.61 ±0.06 



SUMMARY 

We have studied the hysteresis and magnetization pro- 
cesses for N x N arrays (with N = 2 . . . 13) of uniaxial 
ferromagnetic nano dots interacting via the dipole-dipole 
interaction. For an external magnetic field aligned or 
misaligned with one side of the array, the hysteresis loops 
are surprisingly complex. For arrays with odd N the hys- 
teresis loops possess jumps, whereas for even N they do 
not. As N increases, the area A n of the hysteresis loop 
begins to saturate, approaching a non-zero finite value 
determined from a data fit. The area of the hysteresis 
loop scales with N approximately as A^~2 for N odd, 
and approximately as A^~ 2 for even. 
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